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ABSTRACT 

We present numerical relativity results of tidal disruptions of white dwarfs from ultra-close encounters with a 
spinning, intermediate mass black hole. These encounters require a full general relativistic treatment of gravity. 
We show that the disruption process and prompt accretion of the debris strongly depend on the magnitude and 
orientation of the black hole spin. However, the late-time accretion onto the black hole follows the same 
decay, M oc f~ 5//3 , estimated from Newtonian gravity disruption studies. We compute the spectrum of the disk 
formed from the fallback material using a slim disk model. The disk spectrum peaks in the soft X-rays and 
sustains Eddington luminosity for 1-3 yrs after the disruption. For arbitrary black hole spin orientations, the 
disrupted material is scattered away from the orbital plane by relativistic frame dragging, which often leads 
to obscuration of the inner fallback disk by the outflowing debris. The disruption events also yield bursts of 
gravitational radiation with characteristic frequencies of ~ 3.2 Hz and strain amplitudes of ~ 10~ 18 for galactic 
intermediate mass black holes. The optimistic rate of considered ultra-close disruptions is consistent with no 
sources found in ROSAT all-sky survey. The future missions like Wide-Field X-ray Telescope (WFXT) could 
observe dozens of events. 

Subject headings: accretion - black hole physics - gravitational waves - hydrodynamics - relativity - radiation 
mechanisms: general - X-rays: bursts 



1. INTRODUCTION 

Tidal disruptions of stars by black holes (BHs) are fasci- 
nating and violent cosmic events, releasing copious amounts 
of energy in electromagnetic radiation and, in some cases, 
also accompanied by potentially detectable gravitational wave 
emission. The scattered debris from a stellar disruption will 
yield different radiation patterns depending on the orienta- 
tion of the orbit. This variety of radiation patterns will 
provide clues about both the BH and the internal struc- 
ture of the disrupted star. As the debris from a tidal 
disruption showers back into the BH to form an accre- 
tion disk, for instance, it will radiate close to or above 
the Eddington luminosity for some duration, with the spec- 
trum peaking at UV/X-ray frequencies depending on the BH 
mass dReesI [19881 lEvans & Kochaneklll989t ICannizzo et al.l 
Il990t IBogdanovic et al.l 120041: IStrubbe & Ouataerd 120091) . 
In addition, the unbound debris will quickly scatter over 
a large volume, and may be illuminated by the accre- 
tion disk, leadi ng to an optical irradiat i on spectrum domi- 
nated by lines (IBogdanovic et al.l 120041: ISesana et all 120081: 
IStrubbe & Ouataertll2009t iClausen & Eracleousll201 11) . 

When a main sequence star is disrupted, the super- 
Eddington outflow is by itself hot enough to be a sig- 
nificant s ource of optical emission resembling supernova 
radiat ion dKasen & Ramirez-RuTzl 120091: IStrubbe & Ouataerd 
120091) . In these circumstances, instead of outflowing de- 
bris, a steady optic ally thick envelope may reprocess the in - 
ner disk radiation (lLoeb & Ulmerill997l: lUlmer et aDll998l) . 
As a consequence, neither the inner disk nor the irradiated 
debris emission is visible. Another tidal disruption effect 
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arises from the tidal compression perpendicular to the orbital 
plane. The compression leads to the form ation of a strong 
shock that triggers a powerful short outburst ( Kobavashi et alj 
12004 iGuillochon et al.ll2009l; IRosswog et alj|2009l) . In some 
instances, e.g. ultra-close disruptions, the tidal compression 
could be strong enough to produce a thermonuclear ignition 
of the star. In such cases, the detonation of the star would 
be ob served as an under-luminous supernova ( IRosswog et alj 
120091) . 

The variety of proposed radiative signatures often dis- 
agree with each other, calling for precise dynamical model- 
ing of tidal disruptions. E arly hydrodynamic simulations by 
lEvans & Kochanekl (119891) confirmed a flat distribution of de- 
bris mass over its energy range , yielding th e predicted t~ 5 ^ 
law of fallback accretion rate dReesI IT988h . Other numeri- 
cal investigations addressed the tidal compression and related 
shock formation; some of these studie s used smoothed parti- 
cle hydrodynamics (S PH) techniques dKobavashi et alj|2004t 
IRosswog et alj 120091) and ot hers hydrodyna mic grid-based 
codes dGuillochon et al.ll2009l) . The studies by lRosswog et alj 
(120091) included nuc lear reactions with simpl ified networks. 
On the other hand, IBogdanovic et al.l (120041) looked at the 
formation of spherical envelopes lLoeb & Ulmeri ( I 19971) us- 
ing SPH simulations, but the study did not find any con- 
clusive evidence of a spherical structure. At present tidal 
disruption simulations that capture all the relevant physics 
are quite challenging. The simulations require handling, 
in addition to hydrodynamics, effects from general relativ- 
ity, nuclear reactions, and magnetic fields. Thus, many im- 
portant questions are still left unanswered. In particular, 
there are very few simulations that ac count for general rela- 
tivistic effects dLaguna et al.lll993allbllKobayashi et al.ll2004t 
IBogdanovic et all 120041) . The inclusion of magnetic fields 
may lead to jet formation. 

Very likely, the BHs involved in a tidal disruption will have 
spins misaligned with the orbital angular momentum of the 
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incoming star. Thus the disk formed from the tidal debris 
will have angular momentum t hat is also misaligned with the 
BH spin direction dReesll 1988b . For thin disks, the Bardeen- 
Petterson effect will align the BH spin with the angular mo- 
mentu m of the inner region of the disk (Ba rdeen & Pettersonl 
119751) . but for the thick disks ex pected from tidal disruptions , 
the alignment may not happen dPapaloizou & Pringleill983l: 
Fragile et al.ll2007t iStone & Loebll201 11) ." This misalignment 
could have important consequences on the emission channels 
from the disruption, e.g. jet formation. The spin of the BH 
will also produce distinct effects in an ultra-close encounter, 
when the pericentric radius R p is comparable to the gravita- 
tional radius R g = GM^/c 2 of the BH. In this strong grav- 
ity regime, the details of the disruption will depend not only 
on the stellar radius, stellar equation of state, and orbital en- 
ergy, but also on the magnitude and orientation of the spin 
of the BH. In particular, the properties of the accretion will 
depend on the BH spin since its magnitude determines the 
location of the innermost stable circular orbit (ISCO) radius 
Isco, ranging from nsco = 9R g for a counter-rotating maxi- 
mally spinning BH to risco = l^g for a co-rotating maximally 
spinning BH. The value of spin determi nes whether th e star 
following the same orbit gets disrupted ( lKesdenll201 ll) . Fur- 
thermo re, a misaligned rotating BH will drag around n earby 
debris dBardeen et al.ll 19721: [Shapiro & Teukolskvll 19861) . and 
will effectively push this material away from the orbital dis- 
ruption plane. The net effect will be the formation of a shell- 
like disruption debris engulfing the BH, as opposed to the tra- 
ditional S2shap_e_^ebris_oteerved in Newtonian gravity simu- 
lations dEvans & Kochaneklll989l) . 

The present work aims at exploring the exciting regime of 
ultra-close encounters (i.e. R p ~ few R g ) of a white dwarf 
(WD) with an intermediate mass black hole (IMBH). Our 
main goal is to investigate observational signatures due to 
strong gravity effects that may shed light on the presence 
of IMBHs. We consid er a carbon-oxyg en WD with mass 
M wd = 1 Mq and radius dNauenberd 19721) 
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with M ch = 1 .44 M Q the Chandrasekhar mass and R B = 6, 960 
km the Earth mean radius. We fix the mass of the IMBH to 
M hh = 1O 3 M . With these choices, R wd ~ 4R g ~ 6,000km. 
Therefore, the WD and IMBH can be numerically modeled 
with comparable grid resolutions. The other important scale 
in the probl em is the tidal disruption radius R,. For our setup 
dReesll 19881) . 
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gave us ample room to carry out deep penetration encounters. 

Currently, there is a large body of evidence for the existence 
of both solar mass as well as supermassive BHs with masses 
in the range of 10 5 - 10 8 M Q . IMBHs are, however, the miss- 
ing link between stellar mass and supermassive BHs. Today, 
only t entative evidence exists (llrwin et al.H2010l : iFarrell et alj 
l2009HDavis et alj|2011h . On the oth er hand, WDs are thought 
to be abundant in spiral ga laxies dEvans et ai] [19871: 1 eid 
120051) and globular clusters dGerssen et alj|2002l) . Thus the 
identification of distinct signals in both gravitational waves 



(GWs) and the electromagnetic spectrum (Gould 201 1|) from 
a disruption event could potentially guide observations that 
provide evidence for the existence of IMBH as well as in- 
sights into the structure of the WD involved. 

Our study uses the full machinery of numerical relativity 
to solve the Einstein equations of general relativity and hy- 
drodynamics. We do not include nuclear reactions. For the 
ultra-close encounters of our interest, with R p ~ few R g , a 
general relativistic description of gravity is needed. How- 
ever, given the mass ratio M w d/Mbh — 10~ 3 of the bodies in- 
volved and encounters not driven by GW emission, one does 
not need to account for dynamical general relativistic gravity. 
The present study could have been carried out by doing hy- 
drodynamics on the fixed space-time background provided by 
the IMBH. However, this would have implied developing a 
general relativistic hydro code on a fixed background, a code 
similar to the SPH code developed by one of us to investigate 
the ti dal disruption of main sequence stars by supermassive 
BHs dLaguna et al.lll993albl) . We decided instead to take ad- 
vantage of our numerical relativity code Maya. The code has 
demonstrated excellent performance in handling fluid flows in 
the vicinity of BHs in our studies of binary BH mergers in as- 
trophysical envir onments dBode et a l. 2008; Bode et al.ll2.009t 
iBode et ail 1201 11) . The other advantage of using the Maya 
code is the ability to obtain the GW signal directly from the 
simulation, without having to take recourse to the possibly in- 
accurate quadrupole formula or the application of post New- 
tonian (PN) approximations to space-times containing a spin- 
ning BH. 

The paper is organized as follows. In §[2] we describe how 
self-consistent general relativistic initial data is constructed. 
Results from the dynamics of the disruption events are dis- 
cussed in §[3] We estimate electromagnetic transient radiation 
in § fusing a slim disk model to compute the spectrum during 
the fallback phase. We find that the sources shine at Edding- 
ton luminosity ~ 10 41 erg s -1 for about 1-3 yrs, and then fade 
approximately as oc f" 5 / 3 . We also find that there is a 50/50 
chance that the inner disk will be obscured by outflowing de- 
bris for a fully misaligned BH spin. In § [5] we calculate the 
gravitational wave signal produced during the encounters. In 
§[6] we estimate the event rates and discuss the observational 
prospects. We discuss the progress and the limitations of our 
study in §[7] 

2. INITIAL DATA AND CODE TESTS 

When using a fully general relativistic code, such as our 
Maya code, that is based on a 3+1 formulation of the Einstein 
equations, the initial data for the dynamical space-time are 
comprised of the spatial metric 7y and the extrinsic curvature 
Kij of the initial space-time hypersurface. All tensor indices 
are spatial indices and, in this section only, we use units for 
which G=c= 1 . The metric 7,7 characterizes the gravitational 
potentials of the IMBH and WD, and the tensor K,j provides 
the "velocity" of the metric, or the embedding of the space- 
like hypersurface, in the space-timeQ 

Although 7y and K/j are dominated by the IMBH, it is cru- 
cial to account for the contributions from the WD in order to 
correctly include the self-gravity of the WD. The components 
of 7,j and Kjj cannot all be freely specifiable. They need to 

4 For a review on numerical relativity, we recom mend the textbooks by 
lAlcubierrel (2003) and|Baumgarte & Shap^ (M3) 
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satisfy the following elliptic equations 
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Above, is a conformal factor, 5h = diag(l , 1,1) the flat spa- 
tial metric and A its associated Laplace operator. In writing 
Eqs. (0 and (0, we have made the customary assumption 
that the initial physical metric 7y is conformally flat and K,j 
is trace-free. Eqs. (0 and <j4j are, respectively, the conformal 
versions of the so-called Hamiltonian and momentum con- 
straints of general relativity (York 1979). The sources and 
/' are, respectively, the total energy and momentum densities. 

Our approach to solve the Hamiltonian and momentum con- 
straints is to specify the "free" data in Eqs. (0 and (0 from 
the solutions to an isolated, spinning BH, and a boosted WD 
dLoffler et al.ll2006l) . 

Let us consider^ first the solution for a boosted WD. The 
sources and J' in Eqs. (01 and (0 are obtained from 
the stress-energy tensor for a perfect fluid (see chapter 5 in 
iBaumgarte & Shapiroll20 1 Oh . They read 



(9) 
(10) 



p* = p [ ) 



where W = (1 — 7yv ! v-0 -1 ' 2 is the Lorentz factor and v' 
the boost velocity of the WD. In (0 and ([Tol l, p is the 
rest-mass density, P the pressure, and e the specific in- 
ternal energy density. The values for those quantities are 
obtain ed from Tolman-Oppenheime r -Volkoff (TOV ) solu- 
tions dOppenheimer & Volkoffl Il939t iTolmanl fl939l) with a 
poly tropic equation of state P = K p r . We use T = 5/3 appro- 
priate for a non-degenerate gas, which is sufficient to capture 
the dynamics during the inspiral and disruption phases. Dur- 
ing the evolution, we switch to a gamma-law equation of state 
P = pe(T— 1). We denote the solutions to Eqs. (0 and (0 in 
the absence of the BH by wd and Ay d . 

Now we consider the solution to the constraints for a spin- 
ning BH. In the absence of the WD, p* = and J' = 0. In 
this c ase, the momentum co nstraint (0 can be solved analyti- 
cally dBowen & Yorklll980l) . The solution reads 



A*f=-;(e kil S l l k lj + e kjl S l l% 



(11) 



with /' = x'/r and e, ;J t the three-dimensional Levi-Civita sym- 
bol. Above, S' is the spin vector of the BH. 

With the solutions Af? and A^J 1 at hand, it is clear that the 
linear superposition Ay = Ay +Ay d is the solution to the mo- 
mentum constraint (0 in the presence of both a WD and 
a BH._ Next, we solve the Hamiltonian constraint (0 with 
Ay = Ay h +Ay d and given by (0. We use the popular ansatz 



> = 1 + + u. 

2r 



Thus, equation (0 becomes an equation for u. In ([L21 . the 
first term is the flat space solution to (0 ; the second is the BH 
solution; and the third term u includes both the gravitational 
potential of the WD as well as the interaction effects with the 
BH. 

To test our initial data method, and in particular to assess 
the importance of the WD-BH interaction effects, we compute 
initial data for a WD initially at rest; that is, J' = and thus 
Ay d = 0. In addition, we assume a non-spinning BH, i.e. Ay 1 = 
0. The Hamiltonian constraint (0, with the help of ( TTSI i. takes 
the form 

Am = -16tt0- 3 ^. (13) 

The solution to ( TT3l is shown in Figure For color versions 
of the figures in this paper see the electronic edition of the 
Journal. Also in Figure0is the solution m w( j to equation ( [T3l l 
without the BH; that is, with M* = in ([L2"l i. The difference 
between u and k w< j is entirely due to the gravitational interac- 
tion of the WD with the BH. 
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FIG. 1. — Solution u to equation U3t for a WD at rest a distance 15 S w j 
away from a non-spinning IMBH. The solution u provides not only the grav- 
itational potential of the WD but also includes the full non-linear interac- 
tions with the IMBH. For comparison, we plot also i/ w[ ), the solution to equa- 
tion 1 1 3i for a WD in isolation. 

We also tested the ability of our method to produce stable 
isolated WD models and to conserve rest mass. We found 
that, over six dynamical times, the central density of the WD 
did not change more than 1.2% of its initial value and rest 
mass was conserved at the 0.1% level. This is in spite of the 
adaptive mesh refinement infrastructure we use (Carpet) not 
being designed to accurately preserve rest mass conservation 
when interpolating between different refinement levels. 



3. SIMULATION RESULTS 



3.1. 



(12) 



Tidal Disruption Simulation Parameters 

Our study consists of a series of six simulations. As 
mentioned before, the IMBH has a mass Mbh = 10 3 M Q , 
and the WD has a mass M^d = 1 M W and radius ~ 
6,000km (lHamada & SalpeterJfT96l . The WD mass and ra- 
dius correspond to a central density pinMai = 1 .33 x 10 7 g cm -3 , 
T = 5/3, and polytropic constants = 2.41 x 10" 9 (cm 3 g" 1 ) r_1 . 
All the simulations start with the WD in a parabolic orbit, in 
the Newtonian gravity sense, a distance 1.5 R t ~ 60 R g along 
the negative i-axis away from the IMBH. The orbital angular 
momentum of the WD is aligned with the z-axis, as depicted 
in Figure 

Table lists the configurations used in the simulations 
where (3 = R t /R p is the penetration factor for an orbit with 
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FIG. 2. — Configuration setup of our simulations. The WD is located a dis- 
tance 1.5 St ~ 60 R g along the negative .i-axis away from the BH. The orbital 
angular momentum L of the WD is aligned with the £-axis. The orientation of 
the BH spin S is determined by the standard 8 and </> angles in a right handed 
spherical coordinate system. 

pericentric distance R p , and a* = \S\/M 2 is the magnitude of 
the dimensionless spin parameter of the IMBH. The angles 6 
and <j> determine the direction of the BH spin. They are the 
standard angles in a right-handed spherical coordinate system 
as depicted in Figure [2] We use the following convention to 
label the cases we studied. If a case is labeled "BxSy", the 
WD had penetration factor /3 = x. In addition, the value of y 
denotes a non-spinning BH (y=0) or the BH's spin orientation: 
(up, down, in-plane, arbitrary) for y = (u,d,i,a), respectively. 
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P 


a* 


e 




B6S0 
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0.0 


0° 


0° 


B6Su 


6 


0.6 


0° 


0° 


B6Sd 


6 


0.6 


180° 


0° 


B6Si 


6 


0.6 


90° 


0° 


B6Sa 


6 


0.6 


63° 


90° 


B8Sa 


X 


0.6 


63° 


90° 



TABLE 1 

Simulation parameters: = R t /R p is the penetration factor, 

a* = \S\/M 2 IS THE DIMENSIONLESS SPIN PARAMETER OF THE CENTRAL 
BLACK HOLE, 9 AND <fr ARE THE ANGLES BETWEEN S AND THE ORBITAL 
ANGULAR MOMENTUM (£-AXIS) AND THE RADIAL DIRECTION JOINING 
THE IMBH AND THE WD (i-AXIS), RESPECTIVELY. 



In all the simulations, we employ eight levels of mesh re- 
finement with radii 1 .24 R g x 2 l (0 < £ < 7) centered at the BH 
location. In addition to these refinement levels, we surround 
the WD during the pre-disruption phase with five additional 
nested boxes of radii 4.96R g x 2 e (0 < I < 4). The resolu- 
tion on the finest refinement is R g / 19.35, the resolution on 
the level covering the WD is R g /9.615. This mesh refinement 
setup becomes insufficient when the WD begins to be dis- 
rupted. At this point, we turn off the boxes tracking the WD 
and construct a larger mesh with resolution R g /9.615 that in- 
cludes both the IMBH and the highly distorted WD. Once the 
expanding debris has cleared the inner region and the density 
has dropped, we turn off this level again to speed up the sim- 
ulation. 

The mesh refinements affect our ability to conserve rest 
mass throughout the computational domain, in particular if 
the meshes are moving, created, or removed. Figure [3] dis- 
plays the relative gain/loss of rest mass in the computational 
domain, <SM/M wd = (M-M wd )/M wd , where M is the total rela- 
tivistic rest mass and M wd the initial WD rest mass. The com- 
puted value of M takes into account the mass loss through the 
BH horizon and the outer boundary of the computational do- 
main. We find tolerable mass conservation, with violations of 



less than 4% over the course of the simulations. Runs B6Su 
and B6S0 display the worst mass conservation. As we shall 
see, these are the cases in which most of the material of the 
star (> 95%) escapes the BH, forming an expanding cloud 
that reaches the coarsest mesh refinements. Notice also that 
most of the errors in mass conservation accumulate after the 
disruption. Run B6Sd is the only run that loses mass. This 
is the case in which the IMBH accretes almost all of the star 
during the initial passage of the WD. 
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FIG. 3. — Relative gain/loss of rest mass in the computational domain, 
<5M/M w( j = (M-M w( j)/M w( j, where M is the total relativistic rest mass and 
M w< i the initial WD rest mass. The total relativistic rest mass M takes into 
account fluxes into the BH. The only simulation run displaying a net mass 
loss is the case B6Sd in which almost all the material accretes promptly. In 
the other cases, most of the mass gain happens in the short high compression 
phase when the WD is at its periapsis point close to the BH. 



3.2. Tidal compression 

The first stage in a tidal disruption encounter consists of 
the deformations that the incoming star experiences due to 
the tidal forces. In the orbital plane of the star, there are two 
competing tidal deformations. One of them is tidal stretch- 
ing along the radial direction joining the IMBH with the WD. 
The other is tidal compression perpendicular to this radial di- 
rection. In the direction perpendicular to the orbital plane, 
there is also tidal compression. This compression, or squeez- 
ing, is stronger than that in the orbital plane because in the 
latter the radial direction changes rapidly as the star reaches 
periapsis, thus modifying the "contact point" in the star of 
the in-plane compression. We concentrate here on the com- 
pression in the direction perpendicular to the orbital plane. 
This compression c ould, in principle, generate enough heat 
to detonate the star dLuminelll 19851; ILuminet & Marck|[T985l: 
iLuminet & Pichonll 19891) by igniting nuclear reactions at the 
core of the WD dRosswog et alj2009l) . To investigate whether 
a nuclear ignition is likely to happen from tidal compression, 
we plot in Figure [4] the instantaneous maximum temperature 
of the star as a function of the density at that location. To 
a good approximation, this density is also the instantaneous 
maximum density. The instantaneous maximum temperature 
is computed as the mass-density-weighted temperature aver- 
age of the portion of the star whose numerical cells that have 
the largest internal energy and make up for 10% of its mass. 

We calcu late temperatures T from the thermal specific en- 
ergy as (see lLee et all d2005h ) 



m p n mc £th= -k B T(Ni on +N e )n nuc + a BB T 



(14) 
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FIG. 4. — Instantaneous maximum temperature as a function of the density 
at the same location. We compute the density-weighted temperature average 
of the hottest cells in the grid which make up 10% of the total rest mass of 
the WD. The dotted black line separates region s in parameter space wh ere 
the ignition time as estimated from Figure 2 of Dursi & Timmes ( 20061) is 
larger/smaller than the dynamical timescale of the WD, i.e. the region of 
successful explosion from the region where the explosion fails. Densities 
smaller than 10 7 gem -3 are not present in Figure 2 and the estimates for the 
ignition time in this region are based on simple linear extrapolation. In all 
cases investigated, the ignition timescale of the WD is much shorter than 
the dynamical timescale = 1 jyJGp > 0.38 s of the WD. 



with 



^initial r-1 

£th = £-y— j-p , 



(15) 



where the second term in equation ( fl5l l is the "cold" specific 
internal energy in the absence of shocks. N e is the number of 
electrons per nucleon making up the gas, N\ on is the number 
of ions per nucleon, n nllc is the number density of nucleons in 
the gas, m p is the proton (nucleon) mass, ks is the Boltzmann 
constant, and obb here is the radiation constant. For simplicity 
we assume N lon = 1/14, for an equal mixture of oxygen and 
carbon atoms. For N e , which depends on the ionization state 
of the plasma, we use N e = for T < 5 x 10 6 K and N e = 1 /2 
for 5 x 10 6 K < T < 1 x 10 10 K. That is, we ignore any partial 
ionization states of the constituent atoms at low temperatures. 

In all cases investigated, the maximum temperature and 
density displayed in the tracks of Figure |4]reach high enough 
va lues to trigger nuclear burning in the star (cf. Figure 2 
of |D ursi & Timm es (2006)). Interestingly, the temperature- 
de nsity tracks in Figure |4] are qualitatively similar to those 
by iRosswog et al.1 (120091) . who included nuclear networks. 



Run 


/3* 


r max (io 9 K) 


Pmax/ Pinitml 


face 


Ainb 


B6S0 


8.9 


8.6 


7.5 


68% 


19% 


B6Su 


9.44 


6.6 


9.3 


< 1% 


60% 


B6Sd 


11 


12 


3.2 


>99% 


< 0.5% 


B6Si 


7.6 


8.6 


3.3 


65% 


22% 


B6Sa 


9.1 


6.0 


8.1 


2% 


67% 


B8Sa 


10 


9.1 


4.0 


43% 


34% 



TABLE 2 

Actual penetration factor p* measured from the 

SIMULATIONS. MAXIMUM TEMPERATURE Tmax AND COMPRESSION 

p/Pinitial wiTHp 1NIT ,AL = 1.33 x 10 7 gcnr 3 . Fraction / ACC of the 

STAR ACCRETED DURING THE FIRST 2s AFTER DISRUPTION, UNBOUND 
FRACTION / unb AT THE END OF THE SIMULATION AT ~ 6s. 



In Table [21 we list the maximum temperature T max and 
compression Pmax/Pinitiab Here T m . dK and p m3X are the max- 
imum temperature and density attained over the course of 
the evolution. We also list the "actual" penetration factor 



f3* = Rt/R*. Recall that the orbital parameters, and in partic- 
ular R p , used to set the WD in a parabolic orbit were obtained 
using Newtonian gravity. Because of relativistic effects, the 
value of R* will differ from the Newtonian estimate R p , and 
thus j3 In our simulations, we define R* as the distance 
of closest approach of the point within the WD where r max , 
or equivalently p max , is found. Notice that for the B6Su and 
B6Sd cases, r max almost doubles and p max triples when go- 
ing from a spin aligned with the orbital angular momentum 
to one that is anti-aligned. This is because in the anti-aligned 
case (B6Sd), the WD penetrates closer, j3* = 11, than in the 
aligned case (B6Sd), f$* = 9.4. Figure [5] shows the maximum 
temperature r max and maximum density p max as a function 
of the actual penetration factor j3* from Table [2] Because 
of the small range of values covered by j3*, it is not possi- 
ble ^to_wrify_fhe_s£^ng_2}nax oc /3 2 and p max oc /3 3 suggested 
bv lCarter & Luminetl (1 19821) 
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3.3. Outflowing debris and spin 

All cases showed that the mass and orientation of the debris 
after disruption is strongly affected by the BH spin. This is 
because the spin determines the location of ISCO radius, and 
it is also responsible for frame-dragging effects. In particu- 
lar, frame-dragging pulls the material out of the orbital plane, 
thereby changing the orientation of the resulting disk. 

We first discuss three disruptions: disruption by a non- 
spinning IMBH (case B6S0), by an IMBH with spin aligned 
(case B6Su), and by an IMBH with spin anti-aligned (case 
B6Sd) relative to the orbital angular momentum (see Ta- 
ble |2). For the spinning BH cases, the magnitude of the di- 
mensionless spin parameter a* was kept at 0.6. With this 
spin magnitude, nsco = 5,670km for the case B6Su, and 
Hsco = H,600km for t he case B6Sd. In the n on-spinning 
case, nsco = 9,000km dMcClintock et al.ll201 ll) . All three 
cases had penetration factors f3 = R t /R p = 6. From equa- 
tion (0, R t = 4-0 R g = 60,000km, thus a penetration factor of 
6 translates to a pericentric distance R p ~ 10,000km. We 
then expect that the WD in the retrograde case B6Sd will 
pass within the ISCO; in the prograde case B6Su the star will 
mostly stay outside of the ISCO; and in the case B6S0 the 
WD will graze by the ISCO. Because r ISC o > ^ P in the B6Sd 
case, almost all of the star is accreted by the BH soon after the 
disruption. For the cases B6S0 and B6Su, on the other hand, 
most of the stellar debris escapes direct capture, expanding 
away from the hole. The debris cloud has a shape resembling 
a thick circular arch. This arch-like cloud eventually closes 
up and forms a disk around the BH. 
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In all cases, the leading edge of the tidally-disrupted star 
wraps around the BH and crashes into the trailing material 
creating a hot region. For the B6Su, B6Sd and B6S0 cases, 
this hot region is mostly contained within the equatorial plane 
and moves outwards from the central region in the form of a 
spiraling disturbance in the debris disk, eventually appearing 
on the debris surface. Figure [6] shows snapshots in the or- 
bital plane of the density (top row) and temperature (bottom 
row) for runs B6Si and B6Su. The snapshots in the left col- 
umn are for the case B6Si ~ 0.5 seconds after disruption. The 
middle and right columns correspond to case B6Su ~ 0.5 sec- 
onds and ^5.9 seconds after disruption, respectively. In the 
aligned case B6Su, the hot region forms a spiral which moves 
outwards from the point where the leading edge of the WD in- 
tersected with the tail material. At late times this feature has 
been washed out. In run B6Si the hot region is much smaller; 
frame dragging by the BH has pulled material out of the plane 
into a shell surrounding the BH. 

For the cases in which the spin of the BH is not aligned 
with the orbital angular momentum of the binary (i.e. B6Si, 
B6Sa and B8Sa cases), frame dragging is effectively able to 
wrap the debris material around the BH. This effect - which 
is not captured by Newtonian calculations - greatly changes 
the geometry of the debris, blanketing the central region with 
a shell-like cloud of gas. We will address the observational 
consequences of the material surrounding the BH in the next 
section. The effect of the misalignment is evident from Fig- 
ure [6] in the remarkable differences of the spiral patterns ob- 
served between the left and middle snapshots at ~ 0.5 seconds 
after disruption. The left panels are for the B6Si case in which 
the BH spin direction is perpendicular to the orbital angular 
momentum. On the other hand, the middle panels are for the 
B6Su case in which the BH spin is aligned with orbital angu- 
lar momentum. The observed differences are a consequence 
of how the debris responds to the orientation of the BH spin. 

3.4. Prompt and late accretion 

We observe an initial prompt phase of accretion when the 
star swings by the BH. This prompt accretion phase is fol- 
lowed by an intermediate phase in which material initially de- 
flected by the BH slows down and begins to accrete. At late 
times, long after the end of our simulations, the characteristic 
?~ 5//3 fallback behavior takes over. 

Figure [7] shows the flow of matter through the IMBH's 
horizon during the first six seconds after the disruption takes 
place. The instant labeled f t idai is the time at which the star 
enters the tidal radius. The horizontal dashed line labeled 
"noise" denotes the level of accretion due to the atmosphere 
used to model the vacuum regions. We observe that the mag- 
nitude of the prompt accretion rate decreases the more pro- 
grade an orbit is. Both B6S0 and B6Su disrupt at comparable 
distances from the central BH, but the peak accretion rate dif- 
fers by almost two orders of magnit ude. 

As already pointed out in Section [33] the BH spin strongly 
affects what fraction of the star is accreted onto the BH soon 
(~ 2 s) after the star disrupts. Table [2] shows this accreted 
fraction for each case. In the B6Sd case, the star is completely 
accreted. On the other hand, the debris in the B6Su and B6Sa 
cases escapes almost in its entirety during the flyby. Notice 
also from Table [2]that the prompt accretion is also correlated 
with the maximum density /o max the WD is able to reach as a 
consequence of the tidal compression. 

The late-time behavior consists of material raining back 



onto the BH with a characteris tic rate of accretion M oc 
r 5/3 dReeslll988t lPhinney|[T989l) . This rate is maintained 
as long as dMfde^, with e^n being the specific k inetic en- 
erg y, is approximately con stant (lLodato et alJl2009h . We fol- 
low [Ro!!wo£etal] (| 2009l) and compute the fallback time for 
each fluid element at a time long after disruption, based on 
data available at the end of the simulation. At this time, hydro- 
dynamical interactions between the fluid elements are small, 
and each fluid element moves on an almost geodesic orbit. 
Figure [8] shows the results of these calculations. Clearly we 
recover the expected f -5 / 3 behavior for return times longer 
than t > 200s. Notice that earliest times in Figure[8]are com- 
parable to the late time accretion rates found during the simu- 
lation in Figure[7] 

4. ELECTROMAGNETIC SIGNATURES 

Ultra-close (i.e. R p ~ few R g ) tidal disruptions are violent 
events. The conversion of even a small fraction of gravita- 
tional energy into light would result in a powerful electromag- 
netic signature. In this section we will use the hydrodynamics 
results from our simulations to make order of magnitude esti- 
mates of the tidal disruption flare's luminosity, its characteris- 
tic photon energy, and the fallback material's electromagnetic 
signature as it eventually forms a slim accretion disk near the 
BH. 

When the WD descends into the potential well of the BH, 
with the typical dynamical timescale of ~ 1 s, it acquires 
energy in excess of ~ 10 53 erg. Even if a tiny fraction of 
this energy gets converted into radiation on this dynamical 
timescale, it would lead to a huge luminosity which would 
be noticeable at cosmological distances. Furthermore, since 
the WD undergoes a substantial acceleration in the deep grav- 
itational potential of the BH, the streams of matter from 
the disruption will approach speeds V p ~ (GMh/S p )'/ 2 ~ 
c((3 R g /R t y/ 2 ~ 0.4/36C, comparable to the speed of light. 
The streams of tidal debris later collide and transform a por- 
tion of their kinetic energy into thermal energy. Being heated 
by this process in excess of 10 9 K, the dense hot plasma 
quickly produces photons (iKrolik & Piranll20llh . These pho- 
tons, however, are effectively trapped by the debr is, capping 
the lu minosity around the Eddington luminosity (Eddington 



L Edd = 47rCGMhmp =l-4xl0 41 



M h 



10 3 Mr. 



erg s 



(16) 



where m p is the proton mass and oj is the Thomson cross- 
section. Notice that the equilibrium between the gravita- 
tional and spherical radiation fields is reached at a luminos- 
ity L = 2LEdd = 2.8 x 10 41 erg s" 1 due to composition effects. 
The system can by far exceed Ledd in the presence of a colli- 
mated outflow, especially in the first several minutes after the 
disruption. We elaborate on this tidal disruption phenomenon 
in a subsequent paper (Shcherbakov et al. 2012, in prep). In 
the present manuscript we limit ourselves to electromagnetic 
signatures at late times of months to years. 

When the accretion rate falls below the Eddington rate 
MEdd, the accretion becomes radiatively efficient, radiating 
mostly in X-rays with luminosity ~ L Edd . Assuming an ef- 
ficiency of e = 0.1, this luminosity translates into an accretion 
rate 



M E dd = 



0.1c- 



= 2.3 x 10" 5 



( M h 
V10 3 M G 



M Q yr- 



(17) 
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FIG. 6. — Snapshots in the orbital plane of the density (top row, in gem -3 ) and 
left column are for the case B6Si ~ 0.5 seconds after disruption. The middle and 
disruption, respectively. 



temperature (bottom row, in K) for runs B6Si and B6Su. The snapshots in the 
right columns correspond to case B6Su ~ 0.5 seconds and ~ 5.9 seconds after 
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FIG. 7. — Accretion rate during the simulation. ttasj indicates the approxi- 
mate time when the WD enters within the tidal radius R t . The horizontal dot- 
ted black line indicates the mass accretion which we observe initially around 
the black hole. It is a purely numerical artifact from the atmosphere used to 
model the vacuum regions. 



The time elapsed until the accretion rate drops to M^dd can be 
estimated based on the fallback times reported in §[3] We have 
demonstrated that the famous f -5 / 3 fallback law dReesll 19881) 
is also present in our ultra-close encounters. From Figure [8] 
we observe that the accretion rate in all but the B6Sd case 
drops below M^dd at times of ?Edd ~ 1 — 2 yrs. Thus the tidal 
events we are considering can, in principle, shine at Edding- 
ton luminosities for many years. The actual luminosity could 
be substantially sub-Eddington, though, if outflowing debris 
produced while the BH swallows matter at super-Eddington 
rates obscures the accretion flow. In this section we will first 
estimate the disk emission within a slim disk model and then 
consider the added effects of obscuration by the outflowing 
debris. 
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FIG. 8. — Mass accretion rates computed from the estimated fallback times 
of the fluid e lements. The s olid line represents the ?~ 5 / 3 falloff predicted in 
iReel n^88b : |Phimievl {1983). 

4.1. Slim Disk Model 

For our study we employ a slim disk m o del, develope d 
and summarized inlAbramowicz et all dl988l) : lLasotal (119941) ; 
iKawaguchil d2003l) : lSadowskil (120091) . of a simple one-zone 
vertically-integrated dynamical model. We compute modified 
black-body bremsstrahlung emission from the disk. For our 
radiation estimates, we assume that the angular momentum of 
the marginally bound debris accreting at late times, and thus 
the angular momentum of the disk, are aligned with the orbital 
angular momentum of the incoming WD. The disk model will 
then, in principle, be misaligned with the spin of the BH. We 
presume that, despite the misalignment, the slim disk extends 
down to ISCO radius computed for the aligned case. In our 
situation it is challenging to find a more reasonable set of as- 
sumptions for a geometrically thick slim disk. The ultimate 
answers require performing full numerical relativity Simula- 
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tions of the fallback regime, a task requiring modeling dy- 
namical timescales 10 6 times longer than those addressed in 
the present study. 

We take our heavy white dwarf to consist of approximately 
50% carbon and 50% oxygen as suggested by stell ar evo- 
lution computations by iMazzitelli & Danton a (1986). The 
same heavy white dw arf composition wa s assum ed by both 
iRosswog et ail (120091) and lKrolik & Piranl d201 ll) . Thus there 
is 1/14 nuclei and 1 /2 electrons per one nucleon with mass 
~ m p . In terms of nucleon density, the ion and electron den- 
sities are 



— « nuc and 



(18) 



respectively. The corresponding non-relativistic heat capacity 
of neutral and fully-ionized gas per single nucleon of temper- 
ature T are 



Cneutial — ~^kgT 



and 



6 t T 

= jk B T, 



(19) 



respectively. The energy of full ionization is Eqno =< Z 2 > 
!3.6eV = 680eV for the assumed composition, where Z w 7 is 
the charge of the nucleus. A CNO WD is thus fully ionized at 
T > 7ion ~5x 10 6 K. The molecular mass of the fully-ionized 
plasma is p = 1 .75, and its density is p = m p n nuc . 

A fully self-consistent model of a slim disk requires an 
estimate of the disk scale-height H. For a disk in the 
super-Eddington regime, H/r=\ (lAbramowicz et all 119881 : 
lKawaguchill200l IStrubbe & Ouataerj 120091) . In our simu- 
lations, circularization of debris happens within a distance 
3QR g from the BH. We therefore consider a truncated disk 
located between the ISCO radius risco and the maximum ra- 
dius r max = 3QR g . The viscosity parameter is taken to be 
a = 0.1. We solve the equations for the central temperature 
as detailed below, taking into account the effects of viscous 
heating, emission, an d the finite heat capacity of matter. For 
simplicity, we follow IShakura & Sunyaevl (1 19731) and set the 
viscous stress at the ISCO to zero. The energy release per unit 
disk area is then 



F + = 



3r g Mc 2 
87rr 3 



1- 



riscp y /2 



The accretion rate is expressed as 



M = 4-TrrHp c v r , 



(20) 



(21) 



where p c is the one-zone density and v, > is the inflow ve- 
locity. 

We adopt a beta disk model where the viscous stress 
is gi ven by t rr t, = ap s with p g being the gas pres- 
sure dSakimoto & Coronitilll98l l) . A beta disk model is one 
of three options. In the standard [Shakura & Sunyaevl d 1973b 
model, cx p tot = p g + p ra d, whereas in the "mea n field" mod- 
els tff, oc y/pgPtot (see e.g. iDone & Davisll2008l for a review). 
There is no definitive answer yet about which model better 
describes stress in a radiation-dominated accretion disk. The 
beta disk model t r< p = ap g assumption leads to a radial velocity 
of plasma 



ac s = a» 



pm p 



(22) 



where c s is the isothermal sound velocity. The expression ( 1221 
for radial velocity, though not fully self-consistent, helps to 



avoid the infinite increase of t he density near the ISCO inher- 
ent to early slim disk models (lAbramowicz et ai]|1988l) . 

The accretion disk emits predominantly bremsstrahlung ra- 
diation. Since the scattering cross-section in the disk is 
much higher than the absorption cross-section, the emis- 
sion spec trum is described by a modified black body 
spectrum dRvbicki & Li ghtmari 119791: ICzernv & Elvislll987l: 
iKawaguchil 20031) . The photon productio n rate via 
brem sstrahlung is dRvbicki & Lightmanl 1 19791; iKatz et all 
l20T(ih 



2 / m e c A i 

Q=\ -ocfZ n ioa n e a T c\ — rA eff g e ff cm s , 
7T V We 



(23) 



where A e ff ~ g e g ~ 1 , a/ w 0.00729 is the fine structure con- 
stant and a j is the Thomson cross-section. With the high den- 
sities, the production rate Q is high so the radiation quickly 
dominates internal energy and pressure. 

Let us estimate the equilibrium temperature of this 
radiation-dominated flow. The typical nucleon density in the 
simulation varies from n nuc ~ 10 25 cm" 3 near the BH, to a min- 
imum of « m in,nuc ~8 x 10 20 cm" 3 . The typical virial tempera- 
ture is 7i on ~ 10 11 K. As photons are produced, the equilib- 
rium implies 



(24) 



leading to photon temperatures r p h ~ 3 x 10 8 K. The timescale 
for reaching equilibrium can be estimated as 



obbT. 



P h 



1 



•■eq 



ki}T p b Q(T ph ) 



x 10' 



,7/8 



"nuc 



10 _s s, 



(25) 



which is significantly shorter than the dynamical timescale 

?orb ~ 1 S. 

We work in the diffusion approximation of high optical 
depth to find the surface temperature. The corresponding scat- 
tering opacity, r sc 3> 1 , is much larger than the bremsstrahlung 
absorption opacity r a b s . The electron scattering opacity is 



k sc = ajlle w 0.2cm 2 g 1 , 

whereas the absorption opacity is 

(l-e- x )e 6 h 2 Zp c n^: 



(26) 



^abs — 



3cklm e mj,Tj /2 x 3 V k B m e 
4.4 x I0 25 (l-e- x )p c 



(27) 



where x = hv/(k B T) dShapiro & Teukolskvl[T986l) . The total 
emitted power through one side of the disk with surface tem- 
perature T s is then 



a BBC ^4 



K abs 
/"> s 



(28) 



1 sc 



where, in a one-zone model, the central density is taken 
as a proxy for the surface density. Here R denotes 
mean opacities. We use here the Rosseland mean opacity 

dShapiro & Teukolsky|[T986l) . R. dbs = 3.2 x 10 22 Zp f T~ 1/2 . 
The energy flux delivered to the surface is 



fsurf = F- = 



a BB c T; 



(29) 
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To sustain the radiation flow through the surface, this energy 
flux is equal to the radiated flux. Thus from (f28l > and ( f29T > 

we obtain that T s oc T c 16 ^ 9 . In a slim disk, cooling does not 
balance heating, F+ F-. The residual energy is stored within 
the accreting material. The energy balance can then be written 

as 

Q = F + -F-, (30) 
where the residual energy per unit surface area is given as 



emitted. As the accretion rate drops with time, the luminosity 
(the area under the curve) also drops. 



Q = AU+A=-Ha BB T c 4 



1 dT c 1 dpc 
T c dt Or dt 



-IT . (3D 



where dt = —dr/v r . We solve the equations (fJTJ and (1221 for 
the central density p c . The central temperature T c is obtained 
from equation [30] with the help of (|20]l, d29j and (EB- Given 
the central density and temperature, we calculate the surface 
temperature T s from (l28l and ( f29b as well as the correspond- 
ing spectrum. 

Figure [9] shows the dependence of surface temperature T s 
on radius R for the accretion rates inferred from the output of 
the B6Su simulation. Later times correspond to lower accre- 
tion rates. The temperature T s is, in general, higher at higher 
accretion rates, i.e. earlier times. This is because the energy 
release per unit time goes up with accretion rate while the ra- 
diating surface area stays constant. At very low M, our slim 
disk accretion model resem bles the standard thin disk model 
(IShakura & Sunyaevl 19731) so T s -» at the ISCO. In contrast, 
T s at the ISCO saturates at the highest accretion rates. With 
little energy release, F+^> F_ , the energy density is propor- 
tional to the accretion rate Q = F+ oc M. In turn, the scattering 
opacity is also proportional to the accretion rate t sc oc n e oc M. 
The energy density of radiation is proportional to T e 4 . In ad- 
dition there is a factor of M in both the numerator and the 
denominator of the expression ( 1291 for the surface flux. In the 
end, F sm f oc M in the radiatively inefficient regime. 
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FIG. 9. — Surface temperature of the slim disk for the case B6Su at various 
times after the disruption. From top to bottom: at / = 8 months, t = 2 months, 
/ = 3 years, / = 10 years, and t = 30 years, respectively. As time progresses 
the accretion rate decreases from M ~ 4000Me(],j to M ~ 4Me(U- The inner 
disk boundary is chosen to be at nsco = 3.83 Rg. 

The corresponding disk spectra are shown in Figure [10] 
From top to bottom, vL v as a function of energy are shown 
at 2, 8 months, then 3, 10 and 30 years, respectively. Like T s , 
the hard tail of the spectrum also saturates at the highest ac- 
cretion rates. On the other hand, the low-energy tail does not 
saturate at M ~ 4000 Mem- This is because the radiation time 
is still shorter then or comparable to the advection timescale 
at large distances from the BH, where the low-energy tail is 




FIG. 10. — Spectrum for the slim disk formed following the tidal disrup- 
tion at various times. The spectra shown corresponds to the case B6Su as in 
Figure [9] From top to bottom: at t = 2 months, t = 8 months, t = 3 years, 
f = 10 years, and t = 30 years, respectively. The inner boundary is chosen to 
be at nsco = 3.83S g . The spectrum gets harder with time, but the luminosity 
decreases. 

The integrated luminosity is shown in FigureQT]for the case 
B6Su. At late times, i.e. t > 2 years, the luminosity is propor- 
tional to the accretion rate and decreases as f" 5 / 3 . Specifically, 
we estimate a luminosity L ~ 0.05Mc~ for our disk model. 
Notice that the estimated luminosity in a fully general rela- 
tivistic model is L ~ 0.091 Mc 2 for a thin disk without radia- 
tive transfer effects dBardeen et alJI 1972b . At early times t < 
2 years the luminosity saturates at about the Eddington value 
^Edd = 1 .4 x 10 41 erg s" 1 . The luminosity at times t < 2 months 
is not reliable, since the slim disk model may not be an ad- 
equate representation of accretion at rates M > 4000MEdd- 
For instance, when the outflow from the inner disk is present , 
the luminosity may further increase (lOwocki & Gavlevll 19971 : 
IShavivl200UlBegelmarj200HlDotan & Shavivl2011l) . Such a 
computation is, however, beyond the scope of the paper. 
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FIG. 1 1 . — Light curve (solid line) of the slim disk following the tidal dis- 
ruption in the case B6Su. The inner boundary is chosen to be at risco = 
3.83Rg. The luminosity decreases with time, first slowly, then sharply below 
Eddington (horizontal dashed line). The behavior asymptotes to a power-law 
L = ().05Mc 2 oc r 5 / 3 at low luminosity (inclined dashed line), when all ther- 
mal energy is radiated on the way to ISCO. The model cannot be trusted if 
luminosity significantly exceeds Eddington value. 



We consider now the radiative properties of disks from 
the other cases with different penetration parameters j3, spins 
magnitudes a*, and spin orientations. Figure Q~2] shows the 
disk spectra at t = 1 yr after the disruption for all the cases 



10 



Haas, Shcherbakov, Bode, & Laguna 



studied. In all cases, we assumed a slim disk with inner 
boundary at risco and outer boundary at 30 R g . 

The spectra in Figure [12] are different because of two fac- 
tors. The first contributing factor is the accretion rate M, 
which we parameterize by the fallback time. More debris falls 
back onto the BH in the B6Sa case than in the other cases (see 
Figure |8). Therefore the spectrum for the B6Sa case is sys- 
tematically brighter. Similarly, the lower accretion rate for the 
B6Si case makes its spectrum systematically weaker than the 
others. We found that one can achieve coincidence of spectra 
among all a* = 0.6 cases with a time shift. That is, at times 
when the accretion rates M are identical for different simu- 
lations, the spectra are also identical. The spectrum for the 
B6Sa case at about t = 16 months is identical to the spectrum 
for the B6Su case at t = 12 months and to the spectrum from 
the B6Si case at about t = 9 months. Also, more matter is 
swallowed by the BH initially in the B8Sa case when com- 
pared with the B6Sa case, but at later times, less matter that 
will be eventually accreted is scattered. Thus the fall back 
time is lower for the B8Sa case, resulting in a less luminous 
spectrum. The spectra also vary because of changes in the 
ISCO radius. For instance, the ISCO is further from the BH 
in the B6S0 case with spin a* = 0. Since little energy release 
is taking place inside nsco = 6^ g , its spectrum is significantly 
softer. 

4.2. Obscuration 

The edge-on view of a slim disk is self-obscured. The spec- 
trum of a self-obscured disk is not easy to determine. We can 
approximate the spectrum by assuming that the thick outer 
edge at r = 3QR g emits modified black body radiation the way 
the top of the slim disk does. The corresponding edge-on 
spectrum is included in Figure Q~2] as the black dotted line. It 
is softer, since only the material far from the BH contributes. 
The surface area of the thick outer edge is large, though, re- 
sulting in a raised low-energy tail. 
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FIG. 12. — Spectra of our slim disk models at t = 1 yr: B6Su case (top 
solid/black), B6Si case (bottom solid/green), B6Sa case (short-dashed/blue), 
B8Sa (long-dashed/blue), B6S0 case (dot-dashed/red). Differences are due to 
changing spin and fallback time. Dotted line corresponds to an edge-on view 
of a slim disk for b6ap6 simulation. The edge-on disk has a softer spectrum 
and a comparable luminosity. 

Besides exploring the properties of a slim disk, we inves- 
tigated whether the emission from the disk is visible through 
the outflowing debris surrounding the disk. Figure [P3l shows 
the (cos 6,4>) directions in which the debris is scattered in the 
"sky" from the BH point of view. We considered only the 
material outside a sphere of radius 1507? g , corresponding to 
the boundary between the bound and the unbound matter at 



time t ~ 6 s after the disruption. The unbound matter travels 
essentially along straight lines at t w 6 s so Figure Qj] ade- 
quately represents the angular distribution of matter at late 
times. White color in Figure [13] indicates there is no matter 
traveling in that direction whereas all darker colors indicate 
the presence of gas. 

The optical depth is very high at directions that are not de- 
noted by white. The inner disk would be completely obscured 
along those directions. The left panel in Figure [13] shows the 
distribution of the scattered debris for the B6Si case. No- 
tice that although the debris is scattered all over the sky, there 
are some preferential directions manifest as a sinusoid in the 
(cos 0,0) plane. On the other hand, in the case B6Su shown 
on the right panel, the debris is scattered preferentially in the 
disruption plane. Thus the orientation of BH spin determines 
the direction of debris outflow for ultra-close encounters. 

We consider now the situation in which we view a WD- 
BH system with an unknown orientation. Recall that we have 
assumed that the disk at late times is fully aligned with the 
initial disruption plane. As discussed before, the disk can, 
in principle, be viewed both edge-on or face-on, although the 
edge-on disk is self-obscured and has a peculiar spectrum (see 
Figure [T2l. The probability the disk is obscured by the out- 
flowing debris is proportional to the solid angle subtended by 
the debris. This solid angle corresponds to the normalized sur- 
face area on the (cos 0,0) plane (see Figure [T3J. There may 
in principle be four different types of obscuration: i) edge- 
on self-obscured disk also obscured by outflowing debris, ii) 
edge-on self-obscured disk unobscured by outflowing debris, 
iii) face-on disk obscured by outflowing debris, and iv) face- 
on completely unobscured disk. 

When the height of the photosphere equals the disk height 
H = r (inclination 9 = 45°), the probability to encounter a self- 
obscured edge-on disk is p e dge = 71%. The obscuration prob- 
ability by outflowing debris is depicted in Figure [14] for both 
edge-on and face-on disks. The panels show the probability as 
a function of debris location at t = 6 s after the disruption. At 
that time, the boundary between the bound and the unbound 
matter is located at 150^ g . Thus, this obscuration probabil- 
ity is also the obscuration probability p b s for days to years 
later. The obscuration probability ;:> bs|edge of an edge-on disk 
can be inferred from the left panel. Since most debris is scat- 
tered in the disruption plane, this probability is relatively high. 
The values reach p b s | e dge ~ 50%, practically independent of 
spin value, spin orientation or penetration factor. The obscu- 
ration probability of a face-on disk p bs|face can be inferred 
from the right panel. Only the fully misaligned case B6Si 
scatters significant amount of debris perpendicular to the dis- 
ruption plane, which leads to relatively high /? b s |face ~ 40% 
for this simulation. Partially and fully aligned simulations all 
scatter matter within 45° from the disruption plane, and thus 
cannot obscure the face-on disk, so that /5 b s |f ace = 0%. The 
total obscuration probability is 

Pobs = A>bs| face Pface+Pobsl edge (1 ~Pfa.ce)- (32) 

The values are p ^ « 50% for the B6Si case and p ^ w 30% 
for the B6Su, B6Sa, B8Sa and B6S0 cases. The misaligned 
simulation with (3 = 8 yields a marginally higher /? b s |edge = 
52% compared to a similar simulation with (3 = 6 for which 
Pobs|ed g e = 44%. As expected, deeper penetration results in 
wider scattering of debris. However, surprisingly, partially 
aligned and fully aligned simulations with f3 = 6 and a* = 0.6 
yield almost the same obscuration probabilities. 
We have previously stated that if there is debris along the 
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FIG. 13. — Obscuring gas beyond r = \50R g on the (cos0,<^) plane for the fully misaligned spin case B6Si {left panel) and fully aligned spin case B6Su (right 
panel). Here 9 is an inclination angle of the disruption plane and <f> is the azimuthal angle with respect to the BH's point of view (See Fig. [2). Left: The gas in 
this case is scattered to all angles, yet there are substantial holes through which the inner disk could be visible. Right: For this case, the disk will be visible in 
directions away from the disruption plane. 
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FIG. 14. — Obscuration by debris as a function of the debris beyond a radius r ou t from the BH. The obscuration probability is the probability that the inner 
disk is not seen from a random viewing angle. Left: obscuration probability for an edge-on disk with inclination angles above 45° (assuming aligned disk). 
Right: obscuration probability for a face-on disk with inclination angles below 45°. The boundary of bound/unbound matter is around r = 150fi g at t = 6 s 
after disruption. The unbound matter moves approximately radially so the probabilities at r out = 150R g (denoted by a vertical line) correspond to obscuration 
probabilities days to years later. Lines types correspond to case B6Si (upper solid red), B6Su (two lower solid green lines: the darker ones are computed at earlier 
time t = 5.4 s), B6S0 (blue dotted), B6Sd (green dashed), B6Sa (dark dot-dashed), and B8Sa (light dot-dashed). The obscuration probability of an edge-on disk 
is about ?s 50% regardless of spin alignment. Only the B6Sd case, with very little debris, exhibits little edge-on obscuration. In turn, the obscuration probability 
of a face-on disk depends on alignment. The fully misaligned B6Si case is obscured at 50% for nearly face-on disk orientations whereas the partially misaligned 
and aligned cases exhibit no obscuration of a face-on disk. 



line of sight, then the inner disk is obscured. However, the op- 
tical depth of debris decays as the debris expands and eventu- 
ally becomes transparent. Let us estimate the time f trans when 
the optical depth of X-ray absorption equals unity r a b s = 1. 
The actual time will of course depend on the line of sight. 
The debris consists mostly of neutral carbon and oxygen with 
high absorption cross-section er a b s ~2x 10~ 20 cm~ 2 within 
0.5— 10 keV photon energy range dHenkeet al.lll993l) . The 
typical outflow velocity of the debris is a substantial fraction 
of the speed of light v ou t ~ 0.15c. The transparency time is 
then 



/ /unb^wdOabs 

56irm p p obs vl ut 



(33) 



where / un b is the fraction of outbound material. From our 
simulations, we found that / un b ~ 22% (/ un b ~ 60%) for the 



B6Si (B6Su) case. With our obscuration probability estimate 
p obs ~ 50% (p obs ~ 30% ) for the B6Si (B6Su) case, we find 
that ftrans w 2 yr (ftrans ~ 4 yr). The transparency time given by 
equation d33l is an underestimate, since marginally unbound 
matter moves at speeds lower than the average v out . 



5. GRAVITATIONAL WAVE EMISSION 

In addition to producing the electromagnetic signatures 
discussed in previous sections, tidal disruptions events 
are potential sources of GW radiation. A disruption 
event generates a burst of GWs lasting t ~ Rp/V p ~ 
C^/GMbh) 1/2 dKobavashi et alj l200l fRosswog et all 120091) . 
The GW burst produced will have a characteristic frequency 
/ - (GM bh /«p) 1/2 and amplitude h ~ (GM wd R s )/(c 2 R p D), 
with D the distance to the BH. For the case of a WD disrupted 



12 



Haas, Shcherbakov, Bode, & Laguna 



by an IMBH dRosswog et alj|2009h . 

-3/2 

/-0.78/3 3/2 




(34) 



(35) 



These estimates place the GWs from these events at the 
low frequency edge of the advanced LIGO design sensitiv- 
ity range. At this edge, the equivalent strain noise is ~ 
Kr^Hz" 1 / 2 (1LIGO Scientific Collaborationl2009t) . 
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FIG . 15. — Gravitational wave strain as calculated in full general relativity 
(solid red) and from applying the quadrupole formula (dashed blue) for B6Sa 
at a distance of 20 kpc. Depicted are the two polarizations, h+ (top) and h x 
(middle), as well as the amplitude (bottom). 

In Figure[15] we compare the GW strain extracted from the 
spaceti me (solid red) to that obtained using the quadrupole 
formula lRosswog et al.l d20() ( )) applied to the WD's center of 
mass trajectory (dashed blue) for B6Sa. While the qualitative 
behavior of the two polarizations, h+ (top) and h x (middle), 
are similar, the fully nonlinear GW strain differs in both am- 
plitude and phase. 
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FIG. 16. 
of 20 kpc. 



- Gravitational wave strain h from disruption events at a distance 
Depicted are the two polarizations, h+ (solid red line) and h x 



(dashed green line) for cases B6Sa (top) and B8Sa (bottom). 

In Figure [16] we show the GW strain as a function of time 
for disruption events at a distance of 20 kpc. Depicted are 




f(Hz) 

FIG. 17. — Power spectrum of the h+ polarization of GWs for all considered 
cases. We include the advanced LIGO design sensitivity (double-dashed line) 
for comparison. All cases indicate a characteristic frequency of 3.2 Hz. At 
low frequencies, ~ 1 Hz, there is a hint of a dependence of the amplitude on 
the spin orientation. 

the two polarizations of GWs: h+ (solid red line) and h x 
(dashed green line) for cases B6Sa (top) and B8Sa (bottom). 
Time t = denotes approximately the time when the disrup- 
tion takes place. Evident is the burst-like nature of the GWs. 
Notice that both cases exhibit comparable burst duration and 
strain amplitude. This is expected since the only difference 
between B6Sa and B8Sa is their value of (5. Although we 
have stated that case B6Sa has j3 = 6 and B8Sa a /3 = 8, from 
Table [2] we have that the actual penetration factor for both 
cases is (3* ~ 10. Interestingly, the other cases also showed 
comparable burst durations, ~ 4s, suggesting a weak depen- 
dence with the orientation of the spin and penetration factor. 

Figure[l7]shows the power spectrum for all the cases inves- 
tigated where we assume that the events took place at a dis- 
tance of 20 kpc. All cases indicate a characteristic frequency 
of ~ 3.2 Hz. This is about an order of magnitude smaller that 
the estimated frequency using (f34-b with /?* ~ 9, the average 
real penetration factor (see Table 0. This should not be too 
surprising since rough estimates using (f34-b and (l35l l do not 
take into account, among other things, spin effects. Similarly, 
Figure [P71 shows characteristic strain amplitudes of ~ 10~ 18 , 
an order of magnitude higher than the estimate using fl35l l with 
P* ~9. 

At low frequencies, ~ 1 Hz, we note that the B6Sa and 
B6Su cases in particular have higher strain amplitude than the 
rest. Interestingly, these are the same cases which, as seen in 
Fig. [7] exhibit significantly lower prompt accretion in the first 
second after disruption and higher accretion in the ~ 1 - 2 s af- 
ter disruption. The third case with higher prompt accretion in 
the ~ 1 -2 s interval, B8Sa, however, exhibits the lowest low- 
frequency tail. This suggests that relativistic effects create an 
entangled parameter space involving spin and penetration fac- 
tor during the disruption phase. Given that the actual penetra- 
tion factors in all the cases were comparable, 7.6 < (3* < 11, 
further study is required to explore this parameter space. 

Synergistic observations of electromagnetic emission from 
the prompt accretion discussed in Section [4] together with 
GW detections, may provide ways to measuring the spin of 
the BH since in both instances the emission seems to depend, 
although not strongly, on the spin orientation. We are cur- 
rently expanding the parameter space covered by our simu- 
lations to investigate this dependence in more detail. Multi- 
messenger observations will be, of course, only possible if 
the prompt electromagnetic emission is not obscured by the 
debris surrounding the BH. If such obscuration prevents ob- 
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serving prompt accretion, the detection of a GW burst would, 
in that case, still serve as an identifying precursor to the ob- 
servations of the late fallback electromagnetic emission. 

6. RATES AND OBSERVATIONAL PROSPECTS 

6.1. Tidal Disruption Rates 

The disruptions of stars by BHs are rare events. More 
studied are tidal disruptions of stars by SMBHs in galactic 
centers. Thei r rate lO^yr" 1 per galaxy was predicted the- 
oretic ally by dMagorrian & TremaindfT999t I Wang & Merrittl 
120041) and fo und consistent with observations by RO SAT 
dDonlev et aUl2002l), X MM-Newton (lEsquei et al.ll2007l) . and 
GALE X (lGezarill2010l) The volume density of SMBHs is 
about dTundo et al.ll2007l) 



nsMBH 



10" 2 Mpc -3 , 



which leads to tidal disruption rate per unit volume 



Rs 



= 10 2 yr -I Gp C 



(36) 



(37) 



The rate of dis ruption of stars in globu lar cluster (GC) is 
much less certain. IBaumgardt et al.l ( 120041) predicts that, for a 
GC with 10 3 M Q BH, the optimistic disruption rate of stars is 
10~ 7 yr _1 per GC. According to iMcLaughlinl (1 1 9991) . the mean 
GC mass is Mpc = 2.4 • 10 5 M^ , which is larger than clusters 
considered i n IBaumgardt et all (120041) . For a largest cluster 
simulated by IBaumgardt et al.l d2004l) with N = 131, 072 stars 
they find 15% of disrupted stars to be WDs. Therefore the 
optimistic rate of WD-IMBH tidal disruptions is 



D= 1.5-10"V" 



(38) 



per cluster. The numbe r of GCs cha nges dramatically from 
one galaxy to another (lHarrisI [l991h . Howeve r, GCs have 
about constant formation efficiency e c i = 0.0026 ( IMcLaughlinl 
119991) . which is the ratio of GC mass to the initial mass of gas 
in a galaxy. We assume half of baryonic matter in galaxies, 
take the Hubble constant Hq = 70.4km s^Mpc" 1 an d current 
baryonic content = 0.00456 (iKomatsu et al.l20TTl) and find 
the number density of GCs as 



n cc = 0.5 



3// ( 2 Vli,€ c i 
87rG M cc 



■■ 34Mpc" 



(39) 



This density overpredicts the number of GCs in the Milky 
Way, but is consistent with the overpopulated with GCs ellip- 
tical galaxies such as M87. It is not known, what percentage 
of GCs contains an IMBH. For an optimistic estimates, we 
take one 10 3 M Q IMBH per cluster. Then, the disruption rate 
per unit volume is 



Rwd-imbh = 5Q0yi Gpc 



(40) 



This rate is even higher than Rsmbh, but the WD-IMBH dis- 
ruptions are much less luminous, so they are harder to find. 
Also , the rates of tidal dis ruptions in GCs are quite uncer- 
tain (IBaumg ardt et al. 2004]), so even disproving the rate given 
by equation d4Qb does not necessarily disprove the hypothesis 
that every GC has an IMBH. 

Out of all tidal disruptions, we focused on ultra-close ones 
with the ratio of tidal radius to pericenter distance about (3 = 6. 
These constitute a fraction of all disruptions. For a realis- 
tic triaxial potential o f a GC, chaotic feeding would prevail 
dMerritt & Poonll2004l) . Chaotic feeding leads to tidal disrup- 
tion probability proportional to R p for objects coming within 



pericenter distance R p , as a clear consequence of gravitational 
focusing. Therefore, tidal disruptions with f3 > 5 constitute 
1/5 of all events. The correspondent optimistic rate is 



Rwd-imbh,p>5 = 100yr 'Gpc . 
6.2. Observations of X-ray Flares 



(41) 



Considered ultra-close tidal disruptions manifest with a 
year-long X-ray flare at an Eddington luminosity off the cen- 
ter of the galax y. Tidal disruptions by SM BHs rather produce 
an optical flare (IStrubbe & Ouataertl2009l) in the center of the 
galaxy. Let us explore if the optimistic high rates of tidal dis- 
ruptions are consistent with past observations and if the future 
observations have a chance of catching such an event. 

ROSAT survey probed all sky down to Frosat = 4 • 
10~ 13 erg s _1 cm~ 2 flux within 0.1-2.4 keV band (IVoges et al.l 
120001) . The peak of X-ray tidal disruption spectrum at 1 yr 
(see Figure IT2T) is at about 2 keV, so that more than half of 
photons are captured. Taking luminosity L = 4.3 • 10 40 erg s" 1 
of B6Su simulation at t = 1 yr we obtain the maximum dis- 
tance of 

d R0S AT = 25Mpc. (42) 

Therefore, the expected number of observed WD-IMBH dis- 
ruptions is Nrosat = 0.015, which is consistent with no events. 
With Chandra, 100 photons can be obtained over 30 ks from 
the source 200 Mpc away and over 1 Ms from the source 
1 .2 Gpc away, according to the Portable, Interactive Multi- 
Mission Simulator (PIMMS: TMukaill 19931) . However, owing 
to a small field of view of the satellite, the careful choice of 
field is required to observe at least one event. 

In turn, the future mission lik e Wide Field X-ray Telescope 
(WFXT) dConconi et al.ll20loh can detect a substantial num- 
ber of events. It is much more sensitive than ROSAT and 
observes in the same waveband. The flux limits of 3 • 10~ 17 , 
5 • 10 , and 3 • 10 for deep, medium and wide surveys with 
correspondent solid angles 100, 3000, and 20,000 deg 2 yield 
maximum event distances 2.5 Gpc, 600 Mpc, and 245 Mpc. 
This translates into 76, 33, and 16 WD-IMBH disruptions or 
5 times fewer ultra-close disruptions. 

Tidal disruptions can also produce powerful prompt super- 
Eddington X-ray flare with the duration of minutes in addition 
to the Eddington-limited X-ray flux on the timescale of ~ 1 yr. 
The lightcurve, spectrum and detection prospects will be dis- 
cussed in the next paper Shcherbakov et al., (2012, in prep.). 

7. DISCUSSION AND CONCLUSIONS 

Our work addresses disruptions of WDs by IMBHs. This 
is the first fully general-relativistic study that includes both 
strong gravity (the metric of the spinning BH) and dynamical 
gravity effects (GWs). We focused on ultra-close disruptions, 
where the periapsis radius is comparable to the BH gravita- 
tional radius. Our study presents for the first time results on 
the influence of the BH spin on the disruption. 

Our study made some simplifying assumptions. Instead of 
employing a degenerate equation of state for the WD, we used 
an ideal gas law. This is, however, expected to have a very mi- 
nor effect on dynamics. If the entropy is approximately con- 
stant, the evolution of a degeneracy pressure supported star is 
identical to that with an ideal gas pressure. The choice of pres- 
sure prescription does not matter until shocks and turbulence 
take over. Similarly, the absence of explicit radiation pressure 
does not influence the initial stages of the disruption. Cold 
outflowing debris may not have significant radiation pressure 
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support, and only the fallback disk will be radiation pressure 
dominated. Radiation emitted from the inner disk may halt 
the infall of the marginally bound debris. However, the disk 
is self-obscured near the disruption plane where most of de- 
bris is scattered, and the photon mean free path is quite small. 
The photons cannot thus effectively transfer energy from the 
inner flow to the outer flow to unbind the debris. Should radia- 
tion pressure be included, the same internal energy contribute 
a radiation pressure smaller than the gas pressure, and thus the 
total pressure in the medium would be lower. 

Our study does not include nuclear reactions, believed t o 
be important for deep WD disruptions dRosswog et alj|2009h . 
Temperature-density tracks of the matter in our simulations 
as depict ed in Figure [4]are com parable to those found in Fig- 
ure 14 of iRosswog et al.l ( 120091) . The tracks suggest that nu- 
clear burning of carbon and oxygen should occur in the cases 
we studied. The additional energy gained in C and O burn - 
ing into Fe is 1 MeV per nucleon dPursi & Timmesl 120061) . 
This translates to a gain of energy of - 0.001 Mc 2 . On the 
other hand, the descent into the gravitational potential of the 
BH yields energies of ~ 0.1 Mc 2 for a periapsis radius of 
R p ~ 10R S . The dispersion of debris energies after ultra-close 
disruptions is also large. Thus, the extra 0.001 Mc 2 energy 
from the burning will not likely have a dramatic effect on the 
debris. The optical lightcurve may indeed change after an ini- 
tially cold white dwarf is heated by nuclear reactions, but the 
influence on the fallback disk is expected to be small. 

A realistic fallback disk simulation should incorporate mag- 
netic fields. Magnetic fields give rise to large viscosity 
and allow the disk to self-consistently accrete. The as- 
sumed parameter a will be a proxy for viscosity until mag- 
netized simulations can be run for a sufficiently long time. 
The presence of magnetic field s may readily lead to an out- 
flow (Bl andford & Pavndl 19821 ). which may boost luminosity 
far above the Eddington limit or radiation pres sure itself may 
drive an outflow under certain circumstances (lShavivll2001l) . 
The duration of our simulations (~ 7 s) was too limited to 
directly model the fallback. The temporal dependence of ac- 
cretion rate depicted in Figure[8]is a ballistic guess, based on a 
confirmed flat distribution of mass over energy. Any outflow 
or halted inflow at late times may result in deviations of the 
fallback law from M oc t ~ 5 / 3 . 

We make a distinction in the paper between the face-on 
and edge-on disks. In reality, only minor differences between 
these orientations may exist. The disk with extreme accretion 
rate may drive a wind or have a corona, which might com- 
pletely obscure the very inner portions of the disk. In such a 
case, the spectrum would be similar to an edge-on spectrum 
regardless of orientation. The precise quantitative predictions 
of spectra will only be possible with detailed numerical sim- 
ulations of the fallback accretion. 

Our study should be viewed as a step towards realistic tidal 
disruptions of WDs by IMBH. Some of the conclusions may 
hold up despite our limited understanding of matter fallback 
and super-Eddington accretion flows. Our main conclusions 
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are: 

• For a non-spinning BH and a BH with a spin aligned 
or misaligned with the orbital angular momentum, the 
debris after disruption forms a thick accretion disk. 

• For misaligned spins, frame-dragging effects scatter the 
debris around the BH and will often obscure the inner 
region from observation. 

• There is a qualitatively different behavior before and 
after a timescale of ~ 1 yr, when the accretion rate ap- 
proaches Eddington accretion rate. The accretion flow 
luminosity stays around Z,Edd before ~ 1 yr and starts 
dropping as f" 5//3 afterwards. 

• The spectrum peaks at soft X-rays. The spectra are 
similar to thin disk spectra at low accretion rates M <C 

MEdd, 

• Self-obscuration and obscuration by debris are present 
in the system. Self-obscuration leads to softer spectrum 
while obscuration by debris may make a fallback disk 
invisible during most of the active accretion period. 

• The GW signal depends weakly on the orientation of 
the spin. The GW burst will be challenging to be de- 
tected for extragalactic sources. 
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